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We analyzed the magnetoresistivity of a two-dimensional electron system excited 
by microwave radiation in a regime of high intensities and low frequencies. In such 
a regime, recent experiments show that different features appear in the magnetore- 
sistivity response which suggest an anharmonic behavior. These features consist 
mainly in distorted oscillations and new resonance peaks at the subharmonics of the 
cyclotron frequency. We follow the model of microwave-driven electron orbits mo- 
tion which become anharmonic when the ratio of microwave intensity to microwave 
frequency is large enough. 
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Microwave-Induced Resistivity Oscillations (MIRO)^'^ and Zero Resistance States 
(ZRS)^'^ have been continuously attracting much attention in the recent years. From the 
experimental standpoint, different features and improvements are being implemented in 
the corresponding experimental configurations^'^i^'^iSiii^. Thus experimentalist have been 
introducing changes in the orientations, frequencies, polarizations and intensities of the 
different fields involved. Accordingly the results obtained with these configurations are 
getting more and more striking and difficult to explain. In some cases totally unexpected 
outcomes have been obtained^. All these experimental evidences and results establish 
real challenges which can be regarded as crucial tests for the theoretical models currently 
available^»^»^>^»^>^>^. Recently more experimental outcomes have joined the group of 
striking experimental contributions^^i^. In this case MIRO's are specifically studied in 
a regime of high MW intensities and low MW frequencies {w). The remarkable results 
show distorted profiles in the diagonal magnetoresistivity {pxx) oscillations and new reso- 
nance peaks at the subharmonics of the cyclotron frequencies. These effects become more 
important when the ratio MW power to MW frequency is larger and exceed a certain 
threshold. It is also noticeable that the main oscillations are shifted to lower magnetic 
field as the latter ratio increases. 

In this letter we report on a theoretical explanation to this unexpected behavior of pxx 
response obtained in regime of large MW power to w ratio. We follow the model of MW 
driven Larmor orbits^^ under this extreme regime. According to this model, the electronic 
orbit guiding center performs under MW excitation, an oscillating and harmonic motion 
with the same frequency as MW. Thus, any change in the experimental set-up affecting 
the electronic orbit motion, will be reflected in the corresponding pxx response. A very 
direct example is the experimental proposal regarding pxx response to bichromatic MW 
radiation^ and the corresponding theoretical explanation^^. Following this model, if we 
increase the MW power and decrease w, the amplitude of the Larmor orbit oscillation 
becomes larger and larger making eventually the oscillating motion anharmonic. Then 
the electrons in their Larmor orbits become driven anharmonic oscillators. In this letter 
we consider only the case of slightly driven anharmonic behavior. It is well-known that 
a small amplitude driven anharmonic oscillating system has the following characteristics: 
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(1) the forced oscillations contains harmonics that are not present in the driving force; 

(2) subsidiary resonances occur at driving frequencies which are subharmonic of the main 
resonance frequency; and (3) the main oscillations are moved to lower frequencies^i. All 
this features are clearly present in the experimental resultsl^,. A possible microscopic 
mechanism for the anharmonic oscillating motion involves an interplay between a large 
oscillation amplitude and the interaction between electrons in their MW-driven oscillating 
Larmor orbits and the lattice ions. 

Following the MW driven Larmor orbits model, we first obtain the exact expression of 
the electronic wave vector for a two-dimensional electron system (2DES) in a perpendic- 
ular and moderate magnetic field B, and MW radiationr^i^^: 

^(x, y, t) = (j)N[{x-X- xci{t)), {y - yci{t)),t] 



xexp— 
n 



^ , dxd , dyd \ , m*Wc{xcix - yciy) , 



Jp{A^)e'P-' (1) 



oo 

X 

p=— oo 

where (pN are analytical solutions for the Schrodinger equation with a two-dimensional 



(2D) parabolic confinement, known as Fock-Darwin states. X is the center of the orbit 
for the electron spiral motion. Xci{t) and yd{t) are the classical solutions for a driven 2D 
harmonic oscillator. For MW radiation polarized along the current direction (x-direction) 
and a harmonic oscillation regime, the expression for Xd is given hy^^: 

it) = , \ cos wt (2) 

m* a/ [w^ — w'^Y + 7^ 

where 7 is a sample dependent damping parameter which affects dramatically the MW- 
driven electronic orbits movement. Along with this movement there occur interactions 
between electrons and lattice ions yielding acoustic phonons and producing a damping 
effect in the electronic motiort^. is the amplitude of the MW field. L is the classical 
lagrangian, and Jp are Bessel functionsi^i2^i2^. 

In the present regime, x^ and y^ are the solutions for the dynamics of a 2D driven 
classical anharmonic oscillator. Since we do not know the exact nature of the anharmonic 
term in the corresponding potential is impossible to solve analytically the classical equa- 
tion of motion^. Nevertheless we can take an alternative approach if we consider that 
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FIG. 1: Calculated magnetoresistivity pxx as a function of B for different MW intensities. MW 
frequency is 17GHz. Vertical dashed lines correspond to the position of the cyclotron resonance 
{w = Wc) and its subharmnonics. It can be observed very clearly the increasing anharmonicity 
of Pxx as the MW power is getting larger. T=1K. 

the oscillating orbits are only slightly anharmonic. Then, although not harmonic, the 
system can be consider still periodic. As for any periodic function, we can try to express 
Xci through a Fourier series and propose a solution like: 

A °° 

Xci{t) = — + '^^[An cos{nwt) + Bn sm(nwt)] (3) 

n=l 

where Aq, An and are the corresponding Fourier coefficients. 

Now we introduce the scattering suffered by the electrons due to charged impurities 
randomly distributed in the sample. To proceed, following the model described in ref- 
erence [16], firstly we calculate the electron-charged impurity scattering rate 1/r (being 
T the scattering time). Secondly we find the average effective distance advanced by the 
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electron in every scattering jump, that in the case of anharmonicity is given by: 

AX*^^ = AX° + v4o/2 + J2^Ar, cos{nwT) + 5„ sm{nwT)] (4) 

n 

where AX" is the effective distance advanced when there is no MW field present^. The 
longitudinal conductivity axx can be calculated: axx oc / dE ^^"^ (fi — //), being fi and 
ff the corresponding distribution functions for the initial and final states respectively and 
E energy. To obtain pxx we use the relation pxx = 2 — where Oxy — ^ and 
(^xx ^ (^xy Finally we can express pxx as being proportional to a sum of Fourier terms 
which reads: 

Pxx oc cos(riwr) + Bn sin(ntt;r)] (5) 

n 

In order to obtain the Fourier terms in pxx, we have carried out a Fourier synthesis pro- 
cess. This process consists in constructing the pxx form by adding together a fundamental 
frequency (which corresponds to the harmonic case) and overtones of different amplitudes. 
Since at this stage it is impossible to obtain analytical expressions for the Fourier coeffi- 
cients, we have introduced phenomenologically the following ones: 

An = OLn 1 o o — (^nCn and 

Bn = Pn eEo ^ _ ^^(j-^ whcrc a„ and /5„ are anharmonicity terms. Their values 

are getting larger as the anharmonicity increases. 

In Fig.l, we show the calculated pxx, as a function of B for different MW intensities 
and w = VJGHz. Vertical dashed lines are positioned at the cyclotron resonance {w = Wc) 
and its subharmnonics. For the three curves presented different developed expressions of 
pxx have been used. Each one shows the increasing anharmonicity for increasing Eq/w 
ratio. As we said above, they have been obtained through a Fourier synthesis process 
where the anharmonicity coefficients have been introduced phenomenologically, keeping 
the number of Fourier terms as small as possible. Thus for the top panel of Fig. 1: 

pxx oc Ci COS wr -|- O.8C2 COS 2wr + O.2C3 cos3wr + 

Ci [0.4 sin 2wT + 0.2 sin 3wr + 0.2 sin 4:Wt] (6) 
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For the middle panel: 

pxx oc Ci cosiyr + 0.4C2Cos2i/;r + O.lCscosSior + 

Ci[0.2sin2wr + 0.1sin3wr] (7) 

And finally for the bottom panel with the lowest Eq/w ratio we recover the harmonic 



response^: 



Pxx oc Ci cos WT (8) 



These figures illustrate how the p^x profile presents increasing anharmonicity as the 
MW power is also increased. Thus, it can be observed clearly the anharmonicity features: 
distorted profile in the pxx oscillations, new resonance peaks at the subharmonics of 
the cyclotron frequencies and finally it is also remarkable that the main oscillations are 
shifted to lower magnetic fields. All this features corresponds unambiguously to a slightly 
anharmonic behavior, i.e., the system amplitude is not very large yet. However when a 
non-linear system is driven with very large amplitude, new vibrational phenomena appear, 
like vibrations in which the motion only repeat itself after two o more driver periods 
leading the systems finally into a chaotic regime^. This latter case is not consider in this 
letter. 

In Fig. 2 we present calculated pxx versus Wc/w for three different frequencies and 
MW intensities. We can observe, as in the experiment!^, how the anharmonicity effects 
increase with decreasing frequency although the MW power also decreases. This suggests 
that in the evolution to anharmonicity, w plays a more important role than MW power. 
Thus, according to a driven oscillating system, we propose that the ratio which govern 
that evolution be proportional to the corresponding driven amplitude: 

eE, Eo 

^ 7^ (9) 
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FIG. 2: Calculated pxx versus Wc/w for three different frequencies and MW intensities. We can 
observe how the anharmonicity effects increase with decreasing frequencies. T=1K. 
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